NON-WEYL RESONANCE ASYMPTOTICS 
FOR QUANTUM GRAPHS 



E. B. DAVIES AND A. PUSHNITSKI 

Abstract. We consider the resonances of a quantum graph Q that consists of 
a compact part with one or more infinite leads attached to it. We discuss the 
leading term of the asymptotics of the number of resonances of Q in a disc of a 
large radius. We call Q a Weyl graph if the coefficient in front of this leading 
term coincides with the volume of the compact part of Q. We give an explicit 
topological criterion for a graph to be Weyl. In the final section we analyze a 
particular example in some detail to explain how the transition from the Weyl to 
the non-Weyl case occurs. 



1. Introduction 

1.1. Quantum graphs. Let Gq be a finite compact metric graph. That is, Qq has 
finitely many edges and each edge is equipped with coordinates (denoted x) that 
identify this edge with a bounded interval of the real line. We choose some subset of 
vertices of Q , to be called external vertices, and attach one or more copies of [0, oo), 
to be called leads, to each external vertex; the point in a lead is thus identified with 
the relevant external vertex. We call the thus extended graph Q. We assume that Q 
has no "tadpoles", i.e. no edge starts and ends at the same vertex; this can always 
be achieved by introducing additional vertices, if necessary. In order to distinguish 
the edges of Qq from the leads, we will call the former the internal edges of Q. 

In L 2 (Q) we consider the self-adjoint operator H — — with the continuity 
condition and the Kirchhoff boundary condition at each vertex of Q; see Section[2]for 
the precise definitions. The metric graph Q equipped with the self-adjoint operator 
H in L 2 {Q) is called the quantum graph. We refer to the surveys (161 ELY] f° r a general 
exposition of quantum graph theory. 

If the set of leads is non-empty, it is easy to show by standard techniques (see 
e.g. |2H Lemma 1]) that the spectrum of H is [0, oo). The operator H may have 
embedded eigenvalues. 

1.2. Resonances of H. The "classical" definition of resonances is 

Definition 1.1. We will say that k G C, k ^ 0, is a resonance of H (or, by a slight 
abuse of terminology, a resonance of Q) if there exists a resonance eigenf unction 
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/ G Lf oc (Q), f ^ 0, which satisfies the equation 

(l.l) - f "( x ) = k 2 f(x), xeg, 

on each edge and lead of Q, is continuous on Q, satisfies the Kirchhoff's boundary 
condition at each vertex of Q and the radiation condition f(x) = f(0)e lkx on each 
lead of Q. We denote the set of all resonances of H by 1Z. 



Any solution to (1.1) on a lead t = [0, oo) satisfies f(x) = ^ne lkx + 'y' e e~ lkx ; the 
above definition requires that there exists a non-zero solution with all coefficients 
7^ vanishing. It is easy to see that all resonances must satisfy Irak < 0; indeed, if 
ko with Im fco > is a resonance then the corresponding resonance eigenfunction is 
in L 2 (Q), so k 2 is an eigenvalue of H, which is impossible since k 2 , ^ [0, oo). As we 
will only be interested in the asymptotics of the number of resonances in large disks, 
we exclude the case k = from further consideration. In the absence of leads, the 
spectrum of H consists of non-negative eigenvalues and k ^ is a resonance if and 
only if k G K. and k 2 is an eigenvalue of H. 

It is well known (see e.g. [HE]) that the above "classical" definition of a resonance 
coincides with the definition via exterior complex scaling (see [TJ [2U [22] )• m the 
complex scaling approach, the resonances of H are identified with the eigenvalues 
of an auxiliary non-self adjoint operator H(i6), 9 G (0,7r). The algebraic multiplicity 
of a resonance is then defined as the algebraic multiplicity of the corresponding 
eigenvalue of H{i6). We discuss this in more detail in Section [2j where we show that 



the multiplicity is independent of 9. In particular, we show (in Theorem 2.3) that 
any k G K, k ^ 0, is a resonance if and only if k 2 is an eigenvalue of H and in this 
case the corresponding multiplicities coincide. 

We define the resonance counting function by 

N(R) = #{A; : k G TZ, \k\ < R}, R>0, 

with the convention that each resonance is counted with its algebraic multiplicity 
taken into account. Note that the set 1Z of resonances is invariant under the sym- 
metry k — y —k, so this method of counting yields, roughly speaking, twice as many 
resonances as one would obtain if one imposed an additional condition Re (k) > 0. 
In particular, in the absence of leads, N(R) equals twice the number of eigenvalues 
A 7^ of H (counting multiplicities) with A < R 2 . 

1.3. Main result. This paper is concerned with the asymptotics of the resonance 
counting function N(R) as R — > oo. We say that Q is a Weyl graph, if 

(1.2) N(R) = -vol{g )R + o{R), zsR^oc, 

where vol(£? ) is the sum of the lengths of the edges of Q . If there are no leads 
then H has pure point spectrum, resonances are identified with eigenvalues of H 



and Weyl's law (1.2) may be proved by Dirichlet-Neumann bracketing. Thus, every 
compact quantum graph is Weyl in our sense. As we show below, in the presence of 
leads this may not be the case. 
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We call an external vertex v of Q balanced if the number of leads attached to v 
equals the number of internal edges attached to v. If v is not balanced, we call it 
unbalanced. Our main result is 

Theorem 1.2. One has 

(1.3) N(R) = -WR + 0(1), asR^oo, 

7T 

where the coefficient W satisfies < W < vol(£?o)- One has W = vol(£?o) if and 
only if every external vertex of Q is unbalanced. 

This theorem shows, in particular, that as the graph becomes larger and more 
complex the failure of Weyl's law becomes increasingly likely in an obvious sense. 



1.4. Discussion. The simplest example of a graph Q with a balanced external ver- 
tex occurs when exactly one lead I and exactly one internal edge e meet at a vertex. 
In this case, one can merge e and i into a new lead; this will not affect the resonances 
of Q but will reduce vol£?o- This already shows that Q cannot be Weyl. Section [6] 
discusses the second most simple example. 



Our proof of Theorem L2 consists of two steps. The first step is to identify the set 
TZ of resonances with the set of zeros of detA(k), where A(k) is a certain analytic 
matrix-valued function. This identification is straightforward, but it has a subtle 
aspect: this is to show that the algebraic multiplicity of a resonance coincides with 
the order of the zero of det A(k). This is done in Sections [3]-[5] by employing a range 
of rather standard techniques of spectral theory, including a resolvent identity which 
involves the Dirichlet-to-Neumann map. 

The function det A(k) turns out to be an exponential polynomial. By a classical 



result (Theorem 3.2), the asymptotics of the zeros of an exponential polynomial 
can be explicitly expressed in terms of the coefficients of this polynomial. Thus, 
the second step of our proof is a direct and completely elementary analysis of the 
matrix A(k) which allows us to relate the required information about the coefficients 
of the polynomial det A(k) to the question of whether the external vertices of Q are 
balanced. This is done in Section [3j 

Resonance asymptotics of Weyl type have been established for compactly sup- 
ported potentials on the real line, a class of super-exponentially decaying potentials 
on the real line, compactly supported potentials on cylinders and Laplace opera- 
tors on surfaces with finite volume hyperbolic cusps in [23 El S E2] respectively. 
The proofs rely upon theorems about the zeros of certain classes of entire functions. 



Likewise, our analysis uses a simple classical result (Theorem 3.2) about zeros of 
exponential polynomials. 

The situation with resonance asymptotics for potential and obstacle scattering in 
Euclidean space in dimensions greater than one and in hyperbolic space is more com- 
plicated and still not fully understood; the current state of knowledge is described 
in [26], |2]. Here we remark only that generically, the resonance asymptotics in the 
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multi-dimensional case is not given by the Weyl formula. We hope that Theorem L2 
can provide some insight to the mult i- dimensional case. 

Resonances for quantum graphs have been discussed in a recent publication [8]. 
Our paper has very little technical content in common with [8], in spite of their 
common theme. 

1.5. Example. In Section [6] we consider the resonances of a particularly simple 
quantum graph which can be described as a circle with two leads attached to it. 



Theorem 1.2| says that if the leads are attached at different points on the circle, the 



corresponding quantum graph is Weyl, and if they are attached at the same point, 
we have a non-Weyl graph. When the two points where the leads are attached move 
closer to each other and eventually coalesce, one observes the transition from the 
Weyl to the non-Weyl case. We study this transition in much detail. We show that 
as the two external vertices get closer, "half" of the resonances move off to infinity. 
In the course of this analysis, we also obtain bounds on the positions of individual 
resonances for this model. 

The same example was recently considered by Exner and Lipovsky [8] subject to 
general boundary conditions that include the Kirchhoff's boundary condition case 
as a singular limit. Although some of their results are broadly similar to ours, none 
of our theorems may be found in [8]. 



2. Resonances via complex scaling 

Here we introduce the necessary notation, recall the definition of resonances via 
the complex scaling procedure and show that the resonances on the real axis coincide 
with the eigenvalues of H. 

2.1. Notation. Let E mt be the set of all internal edges of Q (i.e. the set of all edges 
of Go) and let E ext be the set of all leads; we also denote E = E mt U E ext . The term 
"edge" without an adjective will refer to any element of E. For e G E mt , we denote 
by p(e) the length of e; i.e. an edge e G E mt is identified with the interval [0,p(e)]. 

Let V be the set of all vertices of G, let V ext be the set of all external vertices, and 
let V mt = V \ V ext ; the elements of V mt will be called internal vertices. The degree 
of a vertex v is denoted by d(v ). The number of leads attached to an external vertex 
v is denoted by q(v); we also set q(v) = for v G V" 111 *. 

If an edge or a lead e is attached to a vertex v, we write v G e. If two vertices u, v 
are connected by one or more edges, we write u ~ v . 

We denote by Goo the graph Q with all the internal edges and vertices removed. 
We let xo an d Xoo be the characteristic functions of Go and Goo- 

Let / : G — > C be a function such that the restriction of / onto every edge is 
continuously differentiable. Then for v G V, we denote by N v f the sum of the 
outgoing derivatives of / at v over all edges attached to v. If v is an external vertex, 
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we denote by A^, nt / (resp. N° xt f) the sum of all outgoing derivatives of / at v over 
all internal edges (resp. leads) attached to v. 

Let C(G) be the class of functions / : Q — > C which are continuous on Q \ V ext 
and such that for each external vertex v the function f(x) approaches a limiting 
value (to be denoted by D^f) as x approaches v along any internal edge and f(x) 
approaches another limiting value (to be denoted by D° xt f) as x approaches v along 
any lead. 

For any finite set A, we denote by \A\ the number of elements of A. We will use 
the identity 



Finally, we use the notation C + = {z E C : Imz > 0}. 

2.2. The operator H(x). The domain of the self-adjoint operator H consists of 
all continuous functions / : Q — > C such that the restriction of / onto any e £ £ 
lies in the Sobolev space W|(e), an d / satisfies the Kirchhoff boundary condition 
N v f = on every vertex v of Q. 

For x G R, let U(x) : L (Q) — > L 2 (Q) be the unitary operator which acts as 
identity on L 2 (Qo) and as a dilation on all leads £ = [0, oo): 



This operator admits an analytic continuation to x G C, which we describe below. 

Definition 2.1. For x G C, the operator H(x) in L 2 {Q) acts according to the 
formula 



The domain of H{x) is defined to be the set of all / : Q — > C which satisfy the 
following conditions: 

(i) the restriction of / onto any e G E lies in the Sobolev space Wf (e); 



(2.1) 





(2.4) 





at every internal vertex v; 



(2.5) 
(2.6) 
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In particular, H(0) is the operator called H so far. For complex x, the operator 
H(x) is in general non-selfadjoint. A standard straightforward computation shows 
that for any x G C the operator H(x) is closed and 

(2.7) #(x)* = #(x). 

2.3. Resonances via complex scaling. The following theorem is standard in the 
method of complex scaling; see [U EH |25l [7] : 

Theorem 2.2. The family of operators H(k), x G C, is analytic in the sense of 
Kato (see e.g. [231 Section XII. 2] J and i/je identity 

(2.8) #(x + x ) = U{kq)H(x)U(-h q ), Vx G C, Vx G R 

holds. The essential spectrum of H(h) coincides with the half-line e~ 2l *[0, oo). Lei 
# G (0,7r); i/ien i/ie sector < argA < 27r — 20, A ^ 0, contains no eigenvalues of 
H(i6), and any A ^ in the sector 2ir — 29 < argA < 2ir is an eigenvalue of H(i9) 
if and only if A = k 2 with k G 1Z. 

For completeness, we give the proof in Section |5j 

As 9 G (0, 7r) increases monotonically, the essential spectrum e~ 2l9 [0, oo) of H(i9) 
rotates clockwise, uncovering more and more eigenvalues A. These eigenvalues are 
identified with the resonances k of H via A = k 2 . If A ^ is an eigenvalue of H(i9), 
9 G (0,7r), 2tt — 29 < argA < 2ir, Kato's theory of analytic perturbations implies 
that the eigenvalue and associated Riesz spectral projection depend analytically on 9. 
Noting ( 2.8[ ) and using analytic continuation it follows that the algebraic multiplicity 



of A is independent of 9. It is easy to see directly that the geometric multiplicity of A 
is also independent of 9. The algebraic (resp. geometric) multiplicity of a resonance 
k is defined as the algebraic (resp. geometric) multiplicity of the corresponding 
eigenvalue A = k 2 of H(i9). 

2.4. Resonances on the real line. The geometric multiplicities of resonances will 
not play any role in our analysis. However, we note that for the Schrodinger op- 
erator on the real line, resonances can have arbitrary large algebraic multiplicity 
j!2j . while their geometric multiplicity is always equal to one. This gives an exam- 
ple of resonances with distinct algebraic and geometric multiplicities. It would be 
interesting to see if one can have distinct algebraic and geometric multiplicities of 
resonances for quantum graphs in the situation we are discussing. We have nothing 
to say about this except for the case of the resonances on the real line: 

Theorem 2.3. (i) If k G M, k ^ 0, is a resonance of H then the algebraic and 
geometric multiplicities of k coincide. 
(ii) Any k G M., k ^ 0, is a resonance of H if and only if k 2 is an eigenvalue of 
H and the multiplicity of the resonance k coincides with the multiplicity of the 
eigenvalue k 2 . 

Proof. 1. Let A > be an eigenvalue of H with the eigenfunction /. If i = [0, oo) 
is a lead, then f(x) = ^£e lkx + 7^e"* fcx , x G £, where k 2 = A. Since / G L 2 (£), we 
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conclude that 7^ 
for all 9 and H(i 

(2.9) 



7^ = and so / = on all leads. It follows that / G Dom H(U 
= Xf. This argument proves that 

dimKer(#(i0) - XI) > dimKer(# - XI). 



2. Let / G Ker(H(i8) - XI), X > 0, 9 G (0,7r). Let us prove that / vanishes 
identically on all leads. Let A = k 2 , k > 0. On every lead, we have 

(2.10) f(x) = f(0)exp(ie i9 kx). 

Consider the difference 



(2.11) u (f) 



Go 



\f(x)\ 2 dx-X [ 

JGo 



\f(x)\ 2 dx 



Go 



\f'(x)\ z dx+ / f"(x)f(x)dx. 



Go 



Integrating by parts, we get 



"(/) 



(NTf)Dtf. 



v£V c 



Using the boundary condition ( |2.5[ ) and formula (2.10), we obtain 

w(/) = i* E |^/| 8 5(«)- 



By the definition ( |27Tip of u(f), we h ave Imu(/) = 0. This yields that \D%*f\ = 
on all external vertices v. By ( |2.10 ), it follows that / vanishes identically on all 
leads. 



3. By combining the previous step of the proof with (2.5) and (2.6) we obtain 
D^f = N™f = 0. It follows that for any / G Ker(#(i0) - XI), X > 0, 9 G (0, vr), 
we have / G Dom if and Hf = Xf. This argument also proves that 

(2.12) dimKer(# - XI) > dim Ker(#(i0) - XI). 

4. It remains to prove that if A > is an eigenvalue of H(i8), 9 G (0, it), then 
its algebraic and geometric multiplicities coincide. Suppose this is not the case. 
Then there exist non-zero elements f,g G Dom H(i9) such that H(i9)g = Xg and 
(H(i9) - XI) f = g. 

By step 2 of the proof, g vanishes on all leads. It follows that on all leads the 
function / satisfies (2.10). Next, since g(x) = —f"(x) — Xf(x) on Q , we have 

(2.13) < / \g{x)\ 2 dx = - 



Gu 



(f"(x) + Xf(x))g(x)dx. 



G,< 



Integrating by parts, we get 



(2.14) - / (f(x) + Xf(x))g(x)dx 



G<, 



f(x)(g"(x) + Xg(x))dx 



Go 



v£V c 



v&V c 



Consider the three terms in the r.h.s. of (2.14). The first term vanishes since 
H(i9)g = Xg. Next, since g = on G^, we have D^ xt g = N^g = for any v G V ext . 



8 



DAVIES AND PUSHNITSKI 



By the boundary conditions (2.5) and (2.6) for g it follows that D^g = N^ nt g = 0. 
Thus, the second and third terms in the r.h.s. of (2.14) also vanish. This contradicts 



3. Proof of Theorem 11.21 

Here we describe the resonances as zeros of det A(k), where A(k) is certain entire 
matrix- valued function. Using this characterisation, we prove our main result. 

3.1. Definition of A(k). Fix k G C + . Let C(k) denote the space of all solutions 
/ G L 2 (Q) to — /" = k 2 f on Q without any boundary conditions. The restriction 
of / G C(k) onto any internal edge e has the form f e (x) = a e e lhx + P e e~ tkx , and 
the restriction of / onto any lead I has the form fg(x) = 7^e 4fcx . Thus, dim C(k) = 
2\E int \ + \E cxt \. 

Let us describe in detail the set of all conditions on / G C{k) required to ensure 
that / is a resonance eigenfunction. If f e denotes the restriction of / to an edge e, 
then we can write the continuity conditions at the vertex v as 

(3.1) fe(v) = C v , Ve 3v, 

where ( v G C is an auxiliary variable. We also have the condition 

(3.2) NJ = 0, v G V. 



Writing down conditions (|3.1[), (3.2) for every vertex v G V, we obtain 



n = J2 d ( v ) + \ v \ = 2 \ Eint \ + i^ cxt i + \ V \ 

vev 

conditions. Our variables are ( v , a e , (3 e , ji, altogether we have 

\V\ + dimC(k) = \V\ + 2\E int \ + \E cxt \ = N 

variables. Let (, a, (3, 7 be the sequences of coordinates ( v , a e , (3 e , 7^ of length \V\, 
\E mt \, \E mt \, \E cxt \ respectively, and let v = (C,a,/3,7) T G C^. We may write the 



constraints (3.1), (3.2) in the form Av = 0, where A is an N x N matrix. Each row 



of A relates to one of the constraints, and each constraint is of the form 
(3.3) y-( + a-a + b-(3 + g--f = 0. 



If the constraint is of the form (3.2 ), then y = and a, 6, g all contain a multiplicative 
factor ik which we eliminate before proceeding. The coefficient a e is 0, ±1, or 
±e lfcp( - e \ and the coefficient b e is 0, ±1, or ±e~ tkp ( e \ The coefficient gi is or 1, and 
the coefficient y v is or — 1. 

We have not specified the order of the rows or columns of A(k). However, the 
object of importance in the sequel is the set of zeros of det A(k), and the choice of 
the order of rows or columns of A(k) will not affect this set. 
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3.2. Example. As an example, let us display the matrix A(k) for a graph which 
consists of two vertices v\ and v 2, two edges e\ and e 2 of length p\ and p 2 and a lead 
attached at V\. In this case we have, denoting Zj = e Pi : 



(3.4) A(k) = 



/ 










z 2 


-zi l 




o\ 










1 


1 


-1 


-1 


1 




-1 

















1 




-1 





1 





1 










-1 








1 





1 










-1 


Zl 





zi l 








V 





-1 





z 2 





Z2 1 





3.3. Resonances as zeros of detA(fc). Although A(k) was defined above for k £ 
C+, we see that all elements of A{k) are entire functions of k £ C. Thus, we will 
consider as an entire matrix-valued function of fc. 

In Sections [4]45] we prove 



Theorem 3.1. Any k^ 7^ is a resonance of H if and only if det A(ko) = 0. In this 
case, the algebraic multiplicity of the resonance k coincides with the order of k as 
a zero of det A(k). 

The first part of this theorem is obvious: by the construction of the matrix A, we 
have det A(k Q ) = iff there exists a non-zero resonance eigenfunction / £ C(k ). The 
part concerning multiplicity is less obvious. Unfortunately, we were not able to find a 
completely elementary proof of this part. The proof we give in Sections [3]-[5] involves 
a standard set of techniques of spectral theory of quantum graphs: a resolvent 
identity involving the Dirichlet-to-Neumann map and a certain trace formula. 



By Theorem |3.1[ the question reduces to counting the total multiplicity of zeros 
of the entire function det A(k) in large discs. As it is clear from the structure of the 
matrix A(k), its determinant is an exponential polynomial, i.e. a linear combination 
of the terms of the type e tcrk , a £ M. Thus, we need to discuss the zeros of exponential 
polynomials. 



3.4. Zeros of exponential polynomials. Exponential polynomials are entire 
functions F(k), k £ C, of the form 

n 

(3.5) F(k) = J2 a r eiark > 

r=l 

where a r , a r £ C are constants. The study of the zeros of such polynomials has 
a long history; see e.g. [TH] and references therein. For more recent literature see 
[20J. Some of these results were rediscovered in [SI El [TTJ, where they were used to 
analyze the spectra of non-self-adjoint systems of ODEs and directed finite graphs. 
The asymptotic distribution of the zeros of F depends heavily on the location of the 
extreme points of the convex hull of the set U" =1 {cr r }. 
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We are only interested in the case in which a r are distinct real numbers. We denote 
ct~ = min{o"i, . . . , a n } and a + = max{o"i, . . . , <r n }. For R > we denote by N(R; F) 
the number of zeros of F (counting the orders) in the disc {k G C : \k\ < R}. The 
following classical statement is from [TjJl Theorem 3]. 



Theorem 3.2. Let F be a function of the form ( |3.5 ), where a r are non-zero complex 



numbers and a r are distinct real numbers. Then there exists a constant K < oo such 
that all the zeros of F lie within a strip of the form {k : |Im (z)\ < K}. The counting 
function N(R; F) satisfies 



N(R; F) = —R + 0(1) 

7T 



as R — > +oo. 



3.5. Estimate for N(R; F). Here we prove the first part of the main Theorem 1.2 



Let F(k) = det A(k). From the structure of A(k) it is clear that F(k) is given by 



(3.5) where (Xy ; (J y £1X6 real coefficients. By Theorem 3.2 it suffices to prove that in 
the representation ( |3.5 ) we have 



(3.6) a + < vol(0 o ), a' > -vol(^ ) 



In order to prove (3.6), let us discuss the entries of A{k) in detail. For simplicity 
of notation we will not draw attention in our equations to the fact that all of the 
matrices below depend on k. 

The matrix A has some constant terms and some terms that are exponential in 
k. The term e %kp ^ can only appear in the column associated with the variable a e 
and the term e~ tkp ^ can only appear in the column associated with the variable 
/3 e . The columns associated with the variables ( and 7 contain only constant terms. 
Since the determinant is formed from the products of entries of A where every 
column contributes one entry to each product, we see that the maximum possible 



value for the coefficient a r in (3.5) is attained when every column corresponding to 
the variable a e contributes the term e tkp ^ and every column corresponding to (3 e 
contributes a constant term to the product. The maximal value of a r thus attained 



will be exactly ^2 eeE int p(e) = vol£ - This proves the first inequality in (3.6). The 



second one is proven in the same way by considering the minimal possible value for 



Of course, the coefficients of the terms e ±jfcvo1 ^ ^ in the representation (3.5) for 



det A may well happen to be zero. Theorem 1.2 will be proven if we show that these 



coefficients do not vanish if and only if every external vertex of Q is unbalanced. 



In what follows, for an exponential polynomial F with the representation (3.5) we 
denote by a ± (F) the coefficient a r of the term e l<Trk , a r = ±vol(^o)- 

3.6. Invariance of resonances with respect to a change of orientation. Be- 
fore proceeding with the proof, we need to discuss a minor technical point. Our 
definition of the matrix A(k) assumes that a certain orientation of all internal edges 
of Q is fixed. Suppose we have changed the parameterisation of an internal edge e 
by reversing its orientation. In other words, suppose that instead of the variable 
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x G [0, p(e)] we decided to use the variable x' = p(e) — x. We claim that this change 
will not affect the zeros of det A(k). 

Indeed, let A'(k) be the matrix corresponding to the new parametrization. The 
matrix A'(k) corresponds to the parametrization of solutions / G C(k) on e by 
f(x) = a' e e lkx + (3' e e~ tkx instead of a e e lkx + (3 e e~ %kx . We have 




and thus det A'(k) = - det A(k). 



3.7. Proof of Theorem \1.2 

ternal vertex v of Q is balanced 
vanishes. 



the balanced case. Assume that a particular ex- 



Below we prove that the coefficient a + (detA) 



Let us re-order the rows and columns of A by reference to the vertex v. We 
assume that q internal edges and q leads are attached to v, q > 2. (The case q = 1 is 
trivial because one may then merge the lead with the edge to which it is connected.) 



Using the observation of Section 3.6, we can choose an orientation of these internal 
edges so that they all end at v (i.e. v is identified with the point p(e) of the intervals 
[0, p(e)]). Let the first 2q rows of A be those relating to the conditions (3.1) for the 
vertex v and let the (2q + l)st row be the one relating to the condition (3.2) for 
the vertex v. The ordering of the remaining rows does not matter. Let the first 
2q columns be related to the variables 71, ... , 7 g , oti, . . . ,a q and let the (2q + l)st 
column be related to the variable ( v ; these variables were all defined in Section |2~lj 
The ordering of the remaining columns does not matter. 

We write A in the block form 



(3.7) 



.4 



B C 
D E 



where B is a (2q + 1) x (2q + 1) matrix. For example, in the case q = 2 we have 



(3. 



where z r 



B 



/ 


1 











-1 \ 







1 








-1 










Zt 





-1 













Z2 


-1 


V 


1 


1 


-zi 


-Z 2 


/ 



The determinant is the sum of the products of entries of A where every column 
contributes one entry to each product. In order for the product to be of the type 
a + e lkml ^°\ each column corresponding to a variable a e must contribute the entry 
g ifcp(e) Thus, the constant entries of the columns corresponding to the variables a e 
are irrelevant to our question and can be replaced by zeros; this will not affect the 
value of a + (detv4). Noticing that the columns of D corresponding to the variables 
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71, . . . , 7g and ( v are all zeros, we conclude that 

(B C 
E 

By a general matrix identity, det A Q = det B det E. Finally, a simple row reduction 



shows that det-B = 0; this is easy to see in the case of (3.8). Thus, the coefficient 



a + (det A) vanishes. By (3.6), it follows that a + < vo\Q , as claimed. 

We note (although this is not needed for our proof) that a~ = — vol Qq both in 
the balanced and in the unbalanced case; this will be clear from the next part of the 
proof. 

3.8. Proof of Theorem 11.2b the unbalanced case. Assume that all external 
vertices are unbalanced. We will prove that 

(3.9) a + = vol(So), o- = - vol(0 o ). 



The proof uses the same reduction as (3.7) in Section 3.7, but the details are some 



what more complicated, since now we have to consider all external vertices. 

We label the external vertices by V\,. . . ,v m where m = |\^ ext |. Let Q r denote 
the graph obtained from Q by adding all the leads of Q that have ends in the set 
{v i, . . . , v r }, so that Q m = Q. Let A r denote the constraint matrix A corresponding 
to the graph Q r and let af = a ± (detA r ). 

By the previous reasoning, the graph Q r is Weyl if and only if a+ ^ and a~ ^ 0. 



Our claim (3.9) follows inductively from the following statements: 



1. The graph Q is Weyl. 

2. The coefficient a~ is non-zero for all r. 

3. For all r, if a^_ 1 ^ then a+ ^ 0. 

Item 1 holds because the operator H on Qq has discrete spectrum and no other 
resonances. The eigenvalues obey the Weyl law by a standard variational argument 
using Dirichlet-Neumann bracketing. 

Let us prove item 3. We reorder the rows and columns of A r with reference to v r 
as in Section 3.7 We assume that p internal edges e\,. . . ,e p and q leads £i,. . . , £ q 



are attached to v r , and q ^ p. The first q + p + 1 columns of A r are those relating 
to the variables 71, . . . , 7 g (associated with £i, . . . , £ q ), ai, . . . ,a p (associated with 
ei,. . . ,e p ), and ( r . The first q + p + 1 rows of A r are those relating to the conditions 



(3.1) and (3.2) for the vertex v r . As in Section 3.7, this allows us to write 



(3.10) A r = g % 

where B r is a (q + p+ 1) x (q + p+ 1) matrix. Writing the matrix A r ^i in the same 
way with reference to the same vertex v r , we obtain 



(3.11) A r _ 



Bj-l Cy_l 

-D T ._i E r 
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where -B r _i is a (p + 1) x (p + 1) matrix. In other words, -B r _i, C r -i, -D r -i are 
the matrices B r , C r , D r with relevant q rows and q columns deleted. The deleted 
columns correspond to the variables ji,- ■ ■ ,7 ? , and the deleted rows correspond to 
the conditions (3.1) associated with the leads £i, . . . , l q . Note that the matrix E r 



is the same in (3.10) and (3.11). 



Next, just as in the argument of Section 3.7, we notice that 



aj~ = a + (det B r det E r ) and ajLi = a + (det B r ^\ det E r ). 
Finally, by a simple row reduction we obtain 

(3.12) det B r = (q — p)z\ . . . z p , 

(3.13) det B r _i = (—p)zi . . . z p , 

where Zj = e %kp( - ej \ It follows that a+ and a^_ 1 differ by a non-zero coefficient 
(p ~ q)/p- This proves Item 3. 

Let us prove Item 2. Here the argument follows that of the proof of Item 3, 
only instead of keeping track of the coefficient of e lkvol ^°" > we need to keep track of 
the coefficient of e~ lkml ^°\ and instead of the variables oti,. . . ,a p we consider the 
variables (3i,. . . , (3 p . Instead of the coefficient (q — p) in (3.12) we get (q +p), which 



never vanishes (even if v r is balanced). This proves our claim. 



4. A RESOLVENT IDENTITY AND ITS CONSEQUENCES 



In order to complete the proof of Theorem 1.2, it remains to provide the proof 



of Theorem |3.1| Theorem 4^2 below provides an explicit formula for the difference 
R x (k) — R^){k) in terms of the Dirichlet-to-Neumann map. This leads immediately to 
the trace formula (4.13), which is the key to our proof of Theorem 3.1 in Section [5j 
The formulae obtained in this section are "complex-scaled" versions of resolvent 
identities well known in the theory of boundary value problems, see e.g. fTO] and 
references therein. 

4.1. Dirichlet-to-Neumann map. Throughout this section, we assume that the 
parameter k G C + is fixed. Let C{k) be as defined in Section 3.1 and let Ai(k) = 
C(k) n C(Q). Each / 6 Ai(k) determines a vector £ 6 C' y ' by restriction to V. 
Conversely, every ( e C |y| arises from a function / G Ai(k); this can be seen by 
comparing dim C(k) with the number of constraints imposed by writing f(v) = ( v , 
v G V. Finally, the assumption k G C + implies that only one function / G M(k) 
corresponds to each set of values £ G C^' (otherwise we would have a complex 
eigenvalue of the operator with Dirichlet boundary conditions on all vertices). This 
shows that we may define the Dirichlet-to-Neumann map A(k) : C' y ' — > C' y ' by 

(A(k)() v = N v f 

where / corresponds to £ as described above and N v was defined in Section 2.1 



This map is a well known tool in the spectral theory of boundary value problems 
and has also been used in quantum graph theory [211 HH]- 
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4.2. The functions (p v and formulae for A. Given v G V, let ip v be the function 
in M.{k) that satisfies 

(p v (u) = 5 UV , Wu,v G V. 

The functions ip v are given by the following explicit expressions. Let v G e, e G -E mt 
and identify e with [0,p] where t> corresponds to the point 0. Then 

. , . . , sin Hp — x) 

4.1 p„x = ^- ^, xG0,p=e- 

sin kp 

In the same way, if e G -E ext and f is identified with the point 0, then 

(4.2) <p v (x)=e ikx , xG[0,oo) = e. 

If the dependence on k needs to be emphasized, we will write (p v (x; k) instead of 
(p v (x). 

Lemma 4.1. If k G C + then the map A(k) is invertible. Its matrix entries are given 
by 

(4.3) A uv = 0, ifu^v,u^v; 

Ek 
— — — -, ifu^v,u~v; 



e6-E" 



(4.5) A vv = ikq{v) — k cot(/cp(e)), /or any u G V; 



egB int 



where q(v) was defined in Section 2.1 



Proof. If A(k)( = 0, then the corresponding function / G M(k) C L 2 {Q) satisfies 
the Kirchhoff's boundary condition at every vertex, which implies that / G Domif 
and Hf = k 2 f. Since Spec(if) = [0, oo) and Imfc > 0, this implies that / = 0. 
Therefore A(k) is invertible. 

By the definition of tp v , we have 

A U v = N u (p v . 



The formulae for the matrix entries are obtained by combining this with (4.1) and 



(4.2). □ 



It follows from Lemma 4.1 that A(k) can be extended to a meromorphic function 



of k G C whose poles are all on the real axis, and that for any u, v G V one has 



(4.6) A uv (k) = A vu (k) and A uv (k) = A uv (-k), k G C. 

In the calculations below the expressions A~l will denote the matrix entries of 
(A(k))-K 
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4.3. The complex-scaled version of ip v . We will need a version of the functions 
<f v pertaining to the "complex-scaled" operator H(x). Let k G C + and x G C be 
such that ke x G C + . Given v G V^,we define (p* by 



(p*(x; k) 



<p v (x;k), tfxeQo; 
<p v (0; k)e^ 2 exp(ike x x), if x G £ = [0, oo), £ G £ ext . 

Clearly, ip* is a solution to the equation H(x)(p* = k 2 tp* on every edge of Q. It is 
also straightforward to see that y?^ G C{Q) and y?^ satisfies the boundary condition 



(]2.5|) on every external vertex. For / G C(G), let us denote 

if f G \/ int 
iVi nt / + e- s */ 2 N° xt f, ifveV 



ext. 



It is straightforward to see that 

(4.7) Kv = NZip*, Wu,veV, 

where the Lh.s. depends on k but not on x. Moreover 



(4.8) <p5(x;k) = <p*(x;-k). 

4.4. The resolvent identity. Let Hp be the self-adjoint operator in L 2 (Q) defined 
by Hpf = —f" with a Dirichlet boundary condition at every vertex of Q. Given 
x G C, we define the "complex-scaled" version of Ho as follows; Hd(x) is the 
operator acting in L 2 {Q) defined by 



(H D (x)f)(x) 



-f"(x), ifxeg , 
-e-^f'ix), lixeG^ 

with a Dirichlet boundary condition at every vertex of Q. Of course, ifo(x) splits 
into an orthogonal sum of operators acting on L 2 (e) for all e G E. We see immedi- 
ately that in addition to its essential spectrum e~ 2x [0, oo), the operator H D (x) has 
a discrete set of positive eigenvalues with finite multiplicities. 

We set 

R D {k) = {H D {x) - k 2 I)-\ R"{k) = (H(x) - k 2 I)-\ 

whenever the inverse operators exist. We denote by R K (k;x,y) (resp. R D (k;x,y)), 
x,y G Q, the integral kernel of the resolvent R K (k) (resp. of R D (k)). 

The fact that Hd(x) and H(x) coincide except for different boundary conditions 
at each of the \V\ vertices indicates that the difference of the two resolvents should 
have rank \ V\. Our next theorem makes this explicit. Formulae of this type are well 
known in the theory of boundary value problems; see e.g. [TO] and references therein. 
In the context of graphs, similar considerations have been used in [131 EH EE l2Tj . 

Theorem 4.2. For any k G C + and any x G C, such that ke* G C +; we have 
(4.9) R*(k;x,y)-RUk;x,y) = - £ A~v(k)ip%(x; k)ip£(y; k), 

«,d£V 

for any x,y G Q . 
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Proof. 1. Let R K (k) be the operator in L 2 {Q) with the integral kernel given by 
R"(k;x,y)=RUk;x,y)- £ A.~*(k)ip*(x; k)ip*(y; k). 

u,v£V 

We need to prove that R K (k) is a bounded operator, that it maps L 2 {Q) into 
DomiJ(x) and that the identities 

(4.10) (H{x) - k 2 I)R K (k) = I 

(4.11) R*{k){H{x) - k 2 I) = I 

hold true. First note that since (p% decays exponentially on all leads, the boundedness 

-k). From 



of R"(k) is obvious. Next, using ( [46) ), ( [48] ) one ob tains R K (k)* = R > 
here and (2.7) by taking adjoints we see that (4.11) is equivalent to 



(#(*) - (-fc) 2 )iT(-A;) = J 



which is (4.10) with —k, x instead of fc, x. We note that k G C + , ke K G C + if and 
only if —k G C+, -fee 57 G C+. Thus, ( |4~TT| follows from flUo] ). 

2. It suffices to prove that for a dense set of elements / G L 2 (Q), the inclusion 
R H (k)f G Domif(x) and the identity 

(4.12) (#(x) - k 2 I)R\k)f = f 

hold true. Let / be from the dense set of all continuous functions compactly sup- 
ported on Q and vanishing near all vertices of Q. Let us check that the function 
g = R K (k)f belongs to DomH(x). It is clear that the restriction of g onto any edge 
e of Q belongs to the Sobolev space W 2 (e). Thus, it suffices to check that g belongs 



to C(Q) and satisfies the boundary conditions (2.5) and (2.6). 

Denote go = Rp(k)f. Since g G Domifu(x), g vanishes on all vertices. There- 
fore go lies in C(Q) and satisfies (2.5) at every external vertex v. As mentioned in 



4.3 



Section 

vertex v. Thus, g also has these properties. 



the functions ip* also belong to C(Q) and satisfy (2.5) at every external 



Our next task is to prove that the boundary condition (2.6) is satisfied for the 



function g. Suppose that / is supported on a single edge, which we identify with 
[0,p]. Then the integral kernel of Rp(k) can be explicitly calculated, which gives 

... [ p sin k(p — x) „ . , , 

0o = / — h L f(x)dx. 

Jo sin kp 

Similarly, if / is supported on a lead [0, oo), then a direct calculation shows that 



,2x 



exp(ike*x) f (x)dx . 



Combining this, we see that for any w G V we have 

NZto = [ f(x)<p*(x)dx. 
Jg 
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Using the last identity and (4.7), for any w E V we get: 



K9= [ f(x)<pZ(x)dx- Kv^v [ f(x)<pZ(x)dx = 0, 



u,v£V 



and so the boundary condition (2.6) is satisfied for g. Thus, g £ Domi/(x), as 
required. 



3. It remains to note that the identity (4.12) follows from the fact that Rp is the 
resolvent of Hr>(>c) and the fact that satisfies the equation H(>t)tp* = k 2 ip* on 
every edge and lead of Q. □ 



4.5. A trace formula. The trace formula (4.13) below results by calculating the 



traces of both sides of (4.9). Since the r.h.s. of (4.9) is a finite rank operator, the 



trace is well defined; the fact that the value of (4.13) does not depend on x can be 



proved by complex scaling, but the direct proof is almost as easy. 



The identity (4.13 ) below can be rephrased by saying that the (modified) perturba- 



tion determinant of the pair of operators H(x), Hn(x) equals det A(k). Statements 
of this nature (for x = 0) are well known in the theory of boundary value problems; 
see e.g. [3] and references therein. The key to our proof of Theorem 3.1 will be (4.13) 
and Lemma 5.1, in which det A(k) and det A(k) are related. 



Theorem 4.3. For any k 6 C + and any x e C, such that ke K G C +; we have 
(4.13) Tr(R K (k) -RT D {k)) ■■ 



jdetA(fe) 
2fcdetA(A;)' 



In particular, the l.h.s. is independent of x. 



Proof. 1. Theorem |4. 2 yields 

(4.14) Tr(iT(A;) - RUk)) = - £ A~^(k)a^ v (k), 

where 
(4.15) 

We next compute the coefficients a uv explicitly. If v ^ u and v ^ u then supp ^ D 



<M= / ^ u {x-k) V :(x-k)dx. 

JQ 



supp <p* — and so a uv = 0. If v ^ u and v ~ u then by (4.1 ) 

rp sinkx smk(p(e)—x^ 



0V 



£ 

e£E mt 



o sinkp(e) sinkp(e) 



-dx 



1 sin kp(e) — kp(e) cos kp(e) 



2k 



u,v(Ee 



(sin kp(e)) 2 
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and finally, 



Y r (£) ( smfcx Y dx+ , v) /°°( e -/2 ( ^ x)) 2 rfx 
f-sJo \smkp(e)J HK J J 



1 ^ ^^(e) — cos A;p(e) sin fcp(e) i 
2k ^ (sin kp{e)) 2 + 2k q[V >' 



2. Noting that a uv depend on k but not on x, a direct calculation using (4.3)-(4.5) 
yields 

2~k^ik^ uv = <J ' av ^^' 

It follows that 

Tr(R x (k) - Rum = Y1 A ^^^ Auv{k) 

as required. □ 

5. Proof of Theorems 13.11 and 12.21 

5.1. Calculation of detA(k). Given k E C, we define 
(5.1) 8(k) = Y[ {ksmkp{e)). 

eS-B int 



Let A{k) be the matrix defined in Section 3.1 



Lemma 5.1. For any k G C +; we have the identity 

2 | £ i»*|.| B tat|_| V | 

(5.2) detA{k) = ± klEintl+lvl 5(k)detA(k), 

where the sign ± depends on the ordering of the rows and columns of the matrix 
A(k). 

Proof. 1. Let us order the rows and the columns of A(k) in such a way that the 
first \V\ rows correspond to the conditions N v (u) = 0, and the first \V\ columns 
correspond to the variables (. Then A(k) can be written in the block form as 



(5.3) A 



M 
-N P 



where is the \V\ x \V\ zero matrix and P is a (2|£ int | + |£ cxt |) x (2|£ int | + |£ cxt |) 
matrix. The elements of N are or 1, the elements of M are 0, ±1, ±e ±thp , and the 



elements of P are 0, ±1, or e ±tkp . For example, the matrix (3.4) is written in this 
form. 
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2. Let us reorder the rows of P in such a way that any two constraints associated 
with the continuity conditions at the two endpoints of the same edge follow one 
another. Let us also reorder the columns of P such that each variable e follows the 



corresponding variable a e . For example, the block P of the matrix (3.4) after such 
reordering will be 



/ 1 



z 1 z x 












1 






1 



Z2 z 2 




o\ 





1 J 



In general, after this reordering, P assumes a block-diagonal structure with blocks 
either of size 2x2 with elements 



gifcp g— ikp 



or of size lxl with the element 1. From here it follows that 

(5.4) det P = ± \\ (2i sin(A;p(e))) = ±{2i)\ Ein %-^5{k). 

e£E int 

In particular, since k £ C+, the matrix P is invertible. 

3. By applying the Schur complement method to ( 5.3[ ) one obtains 

(5.5) det A = detPdet(MP _1 iV). 
Let us prove that 

(5.6) ikMP~ l N = A(Jfc). 

Let £ £ C' y ' and let a = P~ 1 N(. The vector a represents a set of parameters a, (3, 
7. Let / G C(k) be the solution with this set of parameters. The equation Pa = N( 
implies that the solution / is continuous on Q and satisfies f(v) = ( v for any vertex 
v. Next, the coordinates of the vector ikMP~ x N( = ikMa are given by 

ik(Ma) v = N v f. 



This shows that ikMa = A(k)(, as required. 
4. By combining ( 5.4[ )— ( |5.6 ) one obtains 



det A(k) = detP(fc) det(M(A;)p- 1 (A;)A^(A;)) 

= ±(2*) |£int| A;- |i?int| (5(£;) det((iA;)~ 1 A(A;)), 



which yields (5.2) immediately. 



□ 
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5.2. Proof of Theorem 3.1[ 1. Let k G C + and let xo an d Xoo be denned as in 
Section 2.1 Clearly, xo-RdWxo is an orthogonal sum of resolvents of the operators 
—d 2 /dx 2 on the intervals (0,p(e)), e G E mt , with the Dirichlet boundary conditions. 
For each such operator we have that (—d 2 /dx 2 — fc 2 ) -1 is trace class and 



Ti(~d 2 /dx 2 - k 2 )- 1 = ((WP) 2 - k 2 )- 1 

n=0 

E 



cot(A;p) 



A 



{k sin(fcp)) 



2k 2 2k J-^k-Trn/p 2k 2 2k y, J 2k(ksm(kp)) 

Summing over all edges, a direct calculation shows that xo-Rd^Xo is a trace class 
operator and 



(5.7) 



Tr(x R D {k)xo) 



d S(k) 



2k5(k) 



2. Let k G C + , ke K G C + . It is easy to see that the resolvent Rp{k) commutes 
with Xo, and that 

XoRdWxo = XqRdWxo- 

Therefore we have 



(5i 



R x (k) - XooRUk)Xoo = R*{k) - R x D {k) + XoR D (k) X o- 



By combining Theorem 4.3 and (5.8), we obtain 

£detA(A0 



(5.9) Tr(R"(k)- X ooRUk)Xc 



iAk) 



tik 



(5(k) detA(fc)) 



2fcdetA(A;) 2kS(k) 



2k5(k) detA(fc) 



Using Lemma 5.1, we then obtain 

(5.10) Tr(R"(k)-XooR%(k)Xc 
for all k G C+ and ke* G C+. 



|^ int | + \V\ _ jdetAjk) 
2kdetA(k) 



2k 2 



3. The r.h.s. of (5.10) is a single- valued meromorphic function of k G C. Let T x {k) 



be the l.h.s. of (5.10). For each fixed x G C, the function r K (k) is meromorphic 
in C with the cut along the line determined by the condition k 2 G a ess (H(>c)) = 
e~ 2 *[0, oo). In other words, r x is meromorphic and single-valued in each of the two 
half-planes lmke K > and lmke x < 0. By the uniqueness of analytic continuation, 



for each x the identity (5.10) extends to all k such that lmke x > 0. 



4. Let k G TZ with the algebraic multiplicity m(k ) > 1 and let 9 G (0, 7r) with 
-9 < axgk < 0. Then Imkoe 10 > and so the identity (5.10) with x = id holds for 



all k near fc . If 7 is a sufficiently small circle with centre at ko, then the multiplicity 
m(ko) equals the rank, or equivalently the trace, of the Riesz spectral projection 



(5.11) 



P 6 (k ) 



2-ni 



R l \k)2kdk. 



NON-WEYL RESONANCE ASYMPTOTICS FOR QUANTUM GRAPHS 21 

Next, since the operator Hu(i9) restricted to -L 2 (£?oo) has no eigenvalues, the oper- 
ator valued function XooRn(k)Xoo is analytic for \mke ie ^ 0. It follows that 



-^-jxooK{k)xoo2kdk = Q. 



By taking the trace of the difference of the last two equations and using (5.10) we 
obtain 

m(k ) = [Tr(R l9 (k)- X ooRD(k)Xoo)2kdk 



2m 



7 



1 fidetA(k) 



2m J 1 det A(k) ' 
Therefore m(k Q ) equals the order of the zero of detA(k) at k = ko, as required. □ 



5.3. Proof of Theorem 2.2[ This theorem is well known and the proof is presented 



here for the sake of completeness. 

1. First note that by Theorem 4.2 the difference of the resolvents of H(x) and 
Hd(>c) is a finite rank operator. By Weyl's theorem on the invariance of the essential 
spectrum under a relatively compact perturbation we obtain 

a ess (#(x)) = a css (H D (x)) = e- 2 "[0,oo). 

2. The fact that the family H(x) is analytic in the sense of Kato follows again 



from Theorem 4.2, since Hjj(x) is analytic in the sense of Kato and each of the 



functions ip% is analytic in x. 



3. The identity (2.8) can be checked by a direct calculation. 

4. Let k G 1Z and let / be the corresponding eigenf unction. For any 9 G (0, 7r) 
with —9 < axgk < 0, let fe be the function defined formally by fe = U(i9)f. More 
precisely, we set fg = fon Q and 

(5.12) f e {x) = f{0)e w/2 exp{ike ie x) 

for x on any lead I = [0,oo). By the choice of 9, we have lmke ld > and so 
fe £ L 2 (Q). A straightforward inspection shows that fe G Dom H(i9) and H(i9)fe = 
k 2 fe- 

5. Conversely, let A e~ 2ie [0, oo) be an eigenvalue of H{i9) for 9 G (0,7r). Write 
A = k 2 with Im ke l6 > 0. Then, for the corresponding eigenfunction g of H(i9) we 
have g(x) = g(0) exp(ike td x) on any lead of Q. A direct inspection shows that g = fe 



in the same sense as (5.12), where / is a resonance eigenfunction. Thus, k G 71 and 
in particular, Imk < 0. It follows that 2n — 29 < axgk 2 < 2m □ 

6. An example 



Here we consider resonances of a particular simple graph Q(c), where c G [0, 1] is 
a certain geometric parameter. The graph Q{c) was also considered in [HI Section 4], 
but with different boundary conditions at the vertices. The graph Q(c) is Weyl for 
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c < 1 and non-Weyl for c = 1. This section has two goals. The first one is to discuss 
the transition between the Weyl and the non-Weyl cases in order to throw new light 
on the failure of the Weyl law. Our second goal is to obtain rigorous bounds on the 
locations of individual resonances of Q(c), which was not addressed in [S]. 



6.1. Definition of Q(c). Given c G [0, 1), we consider the graph Qo(c) which consists 
of two vertices v\ and v 2 and two edges e\ = [0,pi], pi = (1 — c)n, and e 2 = [0,p 2 ], 
Pi — (1 + c)tt. The vertex v 2 is identified with the point of e\ and with the point 
of e 2 , and the vertex v\ is identified with the point pi of e\ and with the point p 2 
of e 2 . Thus, the graph Qo(c) is simply a circle with the circumference yo\Qq(c) = 2tt 
for all c. We attach a lead t\ at v i and a lead £ 2 at v 2 and denote the thus extended 
graph by Q(c). Geometrically, Q(c) is a circle with two leads attached to it. Finally, 
for c = 1, let Q(c) be the circle of length 2tt with two leads attached at the same 
point. 

We will denote by H(c) the operator — 4-% acting in L 2 (Q(c)) subject to the 
usual continuity and Kirchhoff's boundary conditions at the vertices V\ and v 2 . By 



Theorem 1.2, the graph Q(c) is Weyl if and only if c < 1. At the same time, the 
graph can be regarded as the limit of Q(c) as c — > 1 in an obvious geometric 
sense, so we need to explain what happens to resonances as c — > 1. As we will see, 
roughly speaking, half of the resonances of H(c) move off to infinity as c — > 1. We 
will obtain bounds on the curves along which the resonances move as c increases 
from to 1. 



6.2. The matrix A(k,c) for G(c). Let us display the constraints (3.3) correspond- 
ing to the graph Q(c); the matrix A(k, c) will be built up of the rows corresponding 
to these constraints. We denote Zj = e lkpi ' 2 , j = 1,2. The constraints corresponding 
to the vertex v\ are 

(Ri) aiz\ + dizf - Ci = 

(Ri) ol 2 z\ + f3 2 z 2 2 - Ci = 

(R 3 ) 7i - Ci = 

(-R4) -a x zf + (5\z^ 2 - a 2 z\ + /3 2 Z2 2 + 71 = 0. 

The first three constraints above are the continuity conditions, and the last one is 
the requirement that the sum of the outgoing derivatives vanishes. Similarly, the 
constraints corresponding to the vertex v 2 are 

(R 5 ) ^ + ft - C2 = 

(Re) «2 + 02 ~ (2 = 

(Rr) 72 - C2 = 

(R 8 ) oti - (5i + a 2 - (3 2 + 7 2 = 0. 
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We list these constraints in the order Ri, R 5 , R 2 , R&, R3, R7, R4, Rs, and order the 
variables as a±, f3i, a 2 , P2, 7i, 72, Ci> (2- This leads to the matrix 



A(k,c) 



( 




-r 2 














-1 





\ 




1 


1 

















-1 












z 2 


-2 

Z 2 








-1 
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1 











-1 


















1 





-1 
























1 





-1 










-4 


* 2 - 2 


1 













V 


1 


-1 


1 


-1 





1 








/ 



6.3. Calculation of det A(k : c). The graph Q(c) has a reflection symmetry with 
respect to the midpoints of e 1 and e 2 . This allows to decompose the space C(k) 
into the direct sum of the subspaces corresponding to even and odd functions with 
respect to this symmetry. We use this decomposition to represent the matrix A(k, c) 
in a block-diagonal form where the blocks correspond to the even and odd solutions. 
More precisely, let T\ and T 2 be the matrices 



Ti = 



/ 1 






1 




Vo 



1 





-1 









1 





1 








1 





-1 









1 





1 





V 



c 1 













Z2 













1 
1 






z 






1 





-1 



-1 
1 

-Zl 














1 





1 






-z 2 











1 







\ 



■1 / 






1 

-1 






\ 







1 

-1 / 



A straightforward calculation shows that detTi = detT 2 = 16. Next, let A(k,c) = 
T x A(k, c)T 2 ; the reader is invited to check that the matrix A(k) can be written as 

'A n 

Aodd; 



1 cvcn 





with blocks 



A P 



( 2d 







2C 2 




-1 \ 

-1 

1 -1 



A 



odd 



\ -2iSi -2iS 2 1 J 



( 2iS 1 

2iS 2 


y -2d -2C 2 






1 
1 



-1 \ 

-1 
-1 

J 
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where we have used the notation Cj = cos(kpj/2), Sj = sm(kpj/2), j = 1,2. 
Straightforward calculations of det(v4 even ) and det(Axid) now yield 

Theorem 6.1. For all k G C and all c G [0, 1) one has 

det A(k, c) = AF cvcn (k, c)F odd (k, c) 

where 

F even (k,c) = icos(kcrr) + icos(kn) + 2sin(/c7r), 
F dd(k, c) = i cos(kcTr) — i cos(kTi) — 2 sm(kn). 

We will call the zeros of F even (-,c) (resp. of -F dd(-,c)) the even (resp. odd) 
resonances. It is not difficult to check that the resonance eigenfunctions which 
correspond to even/odd resonances are even/odd with respect to the symmetry 



of the graph G(c). By Theorem 2.3, the real even/odd resonances are actually 
eigenvalues of H(c) and therefore we will call them even/odd eigenvalues. 

Finally, it is not difficult to check that the resonances of H(l) are given, as 
expected, by the zeros of det A(k, 1). In fact, in this case we have F odd (k, 1) = 
— 2sin(fc7r) and 

(6.1) F even (k, 1) = 2ie~ ikn ^0 VA; G C. 

Thus, the resonances of H(l) coincide with the solutions to sin(fc7r) = 0, i.e. they 



are given by k G Z. By Theorem 2.3, these resonances (for k ^ 0) coincide with the 



eigenvalues of H{1) and all of them have the multiplicity one. This shows that for 



c = 1 we have the asymptotics (1.3) with W = n = \ vol^o- 



6.4. Locating the odd resonances. 

Theorem 6.2. (i) For any c G [0,1] , any n G Z and any y > one has 

F od d(n + l~iy,c) ^ 0. 
(ii) For any c G [0, 1] and any k = x — iy with y > \x\/y/3 one has F odd (k, c) ^ 0. 

Proof, (i) By an explicit calculation, 

F A&{n+\-iy,c) = % cos((n + i-iy)7rc) + (-l) n sinh(?/7r) -2(-l) n cosh(?/7r) = A+B, 
where 

\A\ = | cos((n + 1/2 - iy)-Kc)\ 

= | cos((n + l/2)7rc) cosh(y7rc) + i sin((n + l/2)7rc) sinh(?/7rc)| 
< cosh(?/7rc) < cosh(?/7r) 

and 

\B\ = 2cosh(?/7r) — sinh(?/7r) = cosh(?/7r) + e~ yn . 

We deduce that 

\F odd (n + l-iy,c)\ > \B\-\A\ > e"^ > 0. 
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(ii) We start by observing that \F ad(k, c)| > 2A — B where 
A = | sin(fc7r)| 



B = | cos(kn) — cos(/c7rc)| 



kn sin(/c7rs) ds 



If u G R and f > then 

sin(w — iv) = sm(u) cosh(t>) — icos(u) sinh(t>). 

Therefore 

sinh(f) < | sin(u — iv)\ < cosh(f). 
We deduce that A > sinh(?/7r) and 

f 1 \k\ \k\ 

B < \k\ir cosh(yns) ds = — (sinh(?/7r) — smh(yns)) < — sinh(?/7r). 

Jc V V 

These bounds imply that 2 A — B > if 2y > \k\, which yields the theorem imme- 
diately. □ 

It follows that all odd resonances are located in the rectangles 

U° n dd = {x-zy:\x-n\<l 0<y< 2 -^}, n e Z. 

The following statement, in combination with Rouchet's theorem, shows that each 
of the rectangles II° dd contains exactly one odd resonance of algebraic multiplicity 
one for all c G [0, 1] . 

Theorem 6.3. If c = there is a resonance of algebraic multiplicity one at k = 
n — zlog(3)/7r for every odd n G Z and an eigenvalue of multiplicity one at k = n 
for every non-zero even n G Z. There is also a resonance of algebraic multiplicity 
one at k = 0. No other odd resonances or eigenvalues exist if c = 0. 

The proof follows from the explicit formula 

iWM) = ^(e^ + 3)(l-e- i/OT ). 

By the implicit function theorem, we obtain that each of the zeros of F dd(s c) is a 
real analytic function of c G [0, 1] with values in II° dd . The set of all odd resonances 
is therefore the union of a sequence of bounded real analytic curves. 

It is interesting to note that each of these resonance curves intersects the real axis, 
thereby (by Theorem 2.3) giving rise to embedded eigenvalues. This happens at 
rational values of c. More precisely, a direct computation shows that -F dd(&, c) = 
for k G R if and only if 

m — n 

k = m + n, c = tor some m,n EN. 

m + n 



Figure H] plots a typical odd resonance curve as c increases from to 1. The curve 

13 5 
7' 7' 7' 



starts at 7 — zlog(3)/7r when c = and passes through 7 when c = =, 1, I, 1. 
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Figure 1. The odd resonance curve in IT° dd 



6.5. Locating the even resonances. 

Theorem 6.4. (i) For any c G [0,1], any n G Z and any y > one has 

F even (n+ | -iy,c) ^ 0. 

(ii) For any c G [0, 1) and any k = x — iy with y > ,°f , one aas -F dd(&, c) 7^ 0. 

Proof, (i) We have 

F even (n + 1/2 - iy, c) = A - B, 



where A, I? are as in the proof of Theorem 6.2 i). The rest of the proof is the same 
as in Theorem 6.2| (i). 



(ii) For any k = x — iy we have 

-e y7r ' c ' + - > cosh(y7rc) > |cos(x7rc) cosh(y7rc) + i sin(x7rc) cosh(y7rc)| 

(6.2) > |cos(x7rc) cosh(y7rc) + i sin(x7rc) sinh(y7rc)| = |cos(/c7rc)| 

and 

(6.3) |icos(jbr) +2sin(/br)| > - |e l ^| - - |e~ i&7r | = -e yn - -e~ y7r . 

2 2 2 2 
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Now suppose F even (fc,c) = 0; then cos(knc) = —icos(kii) — 2sin(fc7r) and therefore, 



combining (6.2) and (6.3), we obtain 

e yn < e yn]c] + 1 + 3e~ yn . 
If y > log(3)/7r or equivalently e y7T > 3 then 

yn < yvr| C | + 2< ^\c\ , _ 2/7T 

3 

A simple manipulation then yields that y < _u) > anc ^ ^ ne re Q u i r ed result follows. 

□ 

It follows that for c G [0, 1) the even resonances are located in the rectangles 

nr n (c) = {x + ly : \x - n\ < §, < y < ^} . 

Just as in the odd case, the following statement shows that for each n G Z and 
c G [0, 1), the rectangle il^ ven (c) contains exactly one resonance. 

Theorem 6.5. If c = there is an even resonance of the algebraic multiplicity one 
at k = n — zlog(3)/7r for every even n G Z and an even eigenvalue of multiplicity 
one at k = n for every non-zero odd n G Z. There are no other even resonances. 

The proof follows from the explicit formula 

F even (A;,0) = -^(e^-3)(l + e- iA;7r ). 

Just as in the odd case, we obtain that the resonances are given by branches of real 
analytic functions of c G [0,1) with values in II^ ven (c). However, in contrast with 
the odd case, the height of the rectangles II^ ven (c) is not uniformly bounded in c. 
Moreover, we have 

Theorem 6.6. Let n G Z and let k n = k n (c) be the unique solution to F even (k, c) = 
with k n (c) G II^ vcn (c). Then lmk n (c) — > — oo as c — >■ 1. 



Proof. Suppose that the conclusion of the theorem is false. Then there exists a se- 
quence c m — > 1 such that Im k n (c m ) is bounded. By passing to a subsequence we can 
assume that k n (c m ) — > k^ G C as m — > oo. This would imply that -F even (&^°, 1) = 



by the joint continuity of the function F even . This is impossible by (6.1). □ 



Formal calculations and numerical analysis suggest that the rate of divergence of 
Imfc n (c) as c — > 1 is logarithmic. Thus, all even resonances move off to infinity and 
this provides partial explanation for the failure of the Weyl law for c — 1. 

As in the odd case, the even resonance curves intersect the real axis for some 
rational values of k. A direct computation shows that F even (A;, c) = for k G E if 
and only if 

m — n 

k = m + n — 1, c— , for some n, m G N. 

m + n — 1 
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Figure [2] plots a typical even resonance curve as c increases from to 1. The curve 
starts at 4 — i log(3) / tc when c = and passes through 4 when c = |, | and diverges 
to 00 as c — > 1. 




Figure 2. The even resonance curve in n^" 311 
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